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ABSTRACT 

The  structural  survivability  of  shipboard  mast/antenna  systems  subjected  to  underwater 
explosion  can  be  ‘‘designed  in,"  through  the  determination  of  the  structural  dynamics  of 
the  masUantenna  system  This  report  details  the  specialized  application  of  accurate  and 
efficient  analytic  methods  for  the  structural  dynamic  design  analysis  of  shipboard 
mast/antenna  systems  Investigated  herein  are  a  class  of  substructuring  methods,  generally 
referred  to  as  component  mode  synthesis  methods,  which  provide  for  the  rapid  calculation 
of  dynamic  response  of  the  mast 'antenna  structural  system  to  weapons  effects 
Additionally,  the  methods  also  provide  for  the  simulation  of  live  fire  testing  The  methods 
allow  the  individual  antennae  and  the  mast  each  to  be  independently  modeled,  arbitrarily 
combined,  and  the  combined  system  dynamic  response  calculated  to  determine  the 
structural  survivability  of  a  proposed  mast'antenna  configuration  This  rapid  and 
“modular"  component-based  analysis  capability  is  specifically  tailored  for  interactive 
computer-aided  design  analysis  of  shipboard  mast/antenna  systems 
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I.  INTRODUCTION 


This  report  documents  the  development  of  analytic  methods  for  maximizing  the 
combat  survivability  of  shipboard  structural  systems  subjected  to  weapons  effects 
Survivability  will  be  improved  through  the  characterization  of  the  mast/antenna  system 
structural  dynamics  and  the  development  of  specialized  design  analysis  tools  for  the 
prediction  and  minimization  of  dynamic  response  due  to  weapons  effects  The  objective  is 
improved  system  combat  survivability 

Additionally,  the  methods  w  ill  be  developed  in  the  context  of  the  analytic  simulation  of 
live  fire  test  and  evaluation  (LFT&E)  for  shipboard  systems  Those  shipboard  structural 
systems  which  undergo  linear  elastic  dynamic  response  due  to  live  fire  effects  can  be 
evaluated  for  live  fire  survivability  using  the  simulation  methods  to  be  developed,  thereby 
eliminating  the  need  for  actual  LFT&E  for  these  systems  Alternatively,  these  simulations 
w  ill  be  of  benefit  in  the  planning  of  actual  live  fire  test  and  evaluation  (LFT&E)  programs 

The  results  herein  focus  on  shipboard  mast/antenna  structures  Shipboard 
mast/antenna  systems  must  be  designed  to  withstand  moderate  to  severe  shock  loading 
induced  by  underwater  explosion  (I'NDEX)  of  conventional  or  nuclear  type  The 

UNDEX  delivers  devastating  forces  to  the  targets  in  the  form  of  incident  shock  wave 

pressure,  gas  bubble  oscillation,  cavitation  closure  pulses,  and  various  reflection  wave 

effects  These  shock-induced  forces  then  propagate  through  the  ship  to  the  various 

systems,  equipment,  and  top-side  structures  including  the  mast  and  antennae  The 
response  of  the  mast  and  antennae  to  the  UNDEX  shock  wave  is  basically  linear  elastic 
and  vibrational  in  nature  The  mast  and  antennae  tend  to  vibrate  at  their  fundamental 
natural  frequency,  or  at  a  low  range  of  natural  frequencies  The  maximum  amplitude  of 
the  vibration  usually  occurs  after  the  shock  wave  passes  the  ship  The  shock  response 
wave  form  is  remarkably  different  at  various  levels  within  the  ship  In  essence,  the  ship 


2 


acts  as  a  low  pass  structural  filter  which  alters  the  characteristics  of  the  propagating  shock 
wave  from  one  possessing  high  frequency  components  to  one  that  contains  relatively  low 
frequency  components  [1]  Thus,  the  shock  survivability  of  the  mast 'antenna  system, 
which  is  located  top  side,  is  a  vibration  problem  in  which  relatively  low  frequency 
equipment  supr  'rt  excitations  are  observed  The  emphasis  on  design  analysis  relates 
directly  to  the  survival  of  the  mission  critical  systems  on  the  platform  The  ability  of  the 
naval  vessel  to  carry  out  its  mission  after  being  subjected  to  an  UNDEX  threat  depends  on 
the  surv  ivability  of  these  systems,  and  specifically  the  mast/antenna  sy  stem 

Combat  survivability  of  new  systems,  such  as  the  mast/antenna  system  can  be 
“designed  in”  by  accounting  for  the  structural  dynamics  of  the  system  during  the  design 
process  The  methods  developed  herein  focus  on  the  structural  dynamics  of  the 
mast  antenna  systems,  so  that  their  combat  survivability  can  be  directly  addressed  in  the 
design  process  Additionally,  the  methods  will  make  possible  the  improvement  of 
survivability  of  existing  systems  For  example,  survivability  can  be  improved  by 
dy  namically  tuning  and  relocating  antennae  based  on  the  application  of  the  methods  to  be 
described 

A.  BACKGROUND 

The  dy  namic  response  of  a  shipboard  antenna  is  dependent  on  the  dy  namic  interaction 
of  the  antenna  with  the  mast  during  response  to  weapons  effect  Large  dynamic  loads  in 
an  antenna  can  result  if  (a)  the  antenna  is  mounted  on  the  mast  at  a  location  with  large 
accelerations  due  to  weapons  effects,  or  (b)  the  antenna  has  its  natural  frequencies  in  close 
proximity  to  the  excited  natural  frequencies  of  the  mast  In  recent  years,  the  Navy  has  had 
frequent  occurrences  of  shipboard  antenna  systems  failing  structurally  after  being 
subjected  to  shock  due  to  weapons  effects  [2]  In  order  to  design  these  structural  systems 
(i  e  mast  and  antenna)  for  minimum  dynamic  response  and  hence  maximum  survivability  , 
the  structural  dynamic  parameters  which  determine  the  dynamic  response  of  the  system 
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must  be  accurate!)  quantified  The  primary  structural  dynamic  parameters  to  be 
determined  are  the  modal /tetrameters  (i  e  ,  natural  frequencies,  mode  shapes,  modal  mass 
and  damping)  of  the  mast  and  various  individual  antennae  The  modal  parameters  are 
required  to  characterize  the  structural  dynamics  of  each  substructure,  e  g  the  mast  and 
each  antenna,  and  hence  characterize  the  dynamics  of  the  combined  structural  system 
Given  an  accurate  coupled  system  analytic  dynamics  model,  weapons-induced  dynamic 
response  can  then  be  predicted,  and  system  designs  can  be  evaluated  and  optimized  with 
respect  to  survivability  The  coupled  system  analvtic  dynamics  model  can  ser\e  as  the 
basis  for  the  computer  simulation  of  LFT&T 

B.  OVERVIEW  OF  THE  SI  BSTRI  CTI  RE  APPROACH  TO  THE  DESIGN 

ANALYSIS  OF  MAST/ANTENNA  SYSTEMS 

The  methods  described  herein  are  directed  at  the  automated  design  analysis  of 
mast  antenna  systems  The  methods  provide  accurate  estimates  of  the  modal  parameters 
for  a  mast 'antenna  structural  system,  and  therefore  vk  ill  provide  accurate  estimates  of  the 
dynamic  response  due  to  weapons  effects  Genera!))  referred  to  as  “component  mode 
synthesis."  these  substructuring  methods  make  use  of  independent  finite  element  models 
for  the  mast  and  each  antenna  In  ordei  to  allow  a  designer  to  rapidly  assess  for 
survivabilit)  a  large  number  of  candidate  mast  antenna  system  designs,  the  methods  are 
computationally  efficient  as  well  as  accurate  With  respect  to  mast 'antenna  systems,  the 
component  mode  synthesis  process  will  allow  a  designer  to  analytically  “install"  the 
various  antenna  models  into  the  mast  model,  and  rapidly  calculate  coupled  mast 'antenna 
system  UNDEX  dynamic  response  When  incorporated  into  a  computer-aided  design 
environment,  the  complexities  of  the  calculation  will  be  transparent  to  the  designer,  and 
w  ill  allow  the  incorporation  of  self-checks  and  protection  against  user  error  and  misuse 
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The  substructure  approach  to  mast/antenna  structural  dynamic  analysis  can  be  briefly 
outlined  as  follows 

•  A  designer  either  finds  the  dynamic  characteristics  of  the  various  antennae  to  be 
installed  from  a  “catalog”  (database)  of  antenna  modal  parameters,  or  calculates 
individual  antenna  modal  parameters  from  a  finite  element  model  of  the  antenna 
The  modal  parameters  of  the  antenna  constitute  the  antenna  dynamic  model 

•  The  various  antennae  dynamic  models  are  analytically  coupled  with  the  mast 
model,  and  the  dynamic  response  of  the  coupled  mast/antenna  system  due  to 
weapons  effects  is  calculated  If  unacceptable  dynamic  response  levels  are 
calculated,  the  various  antennae  models  can  be  rapidly  repositioned  on  the  mast,  or 
exchanged  with  other  antennae,  and  the  new  dynamic  response  calculated 

This  scheme  has  several  significant  advantages  for  the  automated  design  analysis  of 
mast  antenna  systems  The  primary  advantages  include 

•  The  ability  of  these  methods  to  treat  the  mast  and  antennae  as  “substructures,”  and 
arbitrarily  and  repeatedl}  combine  them  for  the  rapid  calculation  of  dynamic 
response  will  make  possible  the  evaluation  of  a  greater  number  of  mast/antenna 
configurations,  and  hence  will  greatly  facilitate  the  determination  of  an  optimal 
configuration  with  respect  to  combat  survivability. 

•  The  various  masts  and  antennae  are  fabricated  by  various  independent  contractors 
The  component  mode  synthesis  method  allow  s  the  separate  modeling  of  the  mast 
and  antennae,  and  therefore  naturally  preserves  the  independence  of  the 
contractors 

•  The  formulations  to  be  described  are  modal,  and  therefore  can  function  equally 
well  with  analytically  derived  modal  parameters,  or  with  modal  parameters 
identified  in  a  vibration  test 
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The  analytic  methods  for  the  generation  of  the  coupled  mast/antenna  model  are  the 
focus  of  this  w  ork  To  be  ev  aluated  in  this  report  are  several  component  mode  methods 
for  substructure  synthesis  the  Craig-Bampton  method  and  two  residua!  flexibility 
formulations  The  methods  are  specialized  for  the  mast/antenna  analysis  problem,  and 
their  relative  merits  compared  in  the  context  of  combat  survivability  The  methods  are 
based  on  the  modal  representation  of  components,  that  is,  rather  than  representing  a 
structure  using  the  mass  and  stiffness  matrices  generated  in  a  finite  element  model,  these 
methods  employ  various  classes  of  “mode  shapes”  to  represent  the  substructures  or 
components  For  example,  the  familiar  normal  modes  of  vibration  are  one  class  of  mode 
shape  used 

The  computational  efficiency  of  these  methods,  which  is  critical  to  their  effectiveness 
in  a  computer  aided  design  environment,  comes  from  their  ability  to  accurately  describe  a 
component  with  a  minimum  number  of  mode  shapes  The  sections  of  this  report  w  hich 
follow  will  describe  the  above  mentioned  synthesis  formulations,  and  demonstrate  their 
relative  accuracy  and  efficiency  in  the  calculation  of  the  dynamic  response  of  a  small  yet 
representativ  e  mast'antenna  model,  subjected  to  a  variety  of  applied  harmonic  forces  as 
well  as  deck  accelerations  and  displacements  The  model  used,  which  includes  a  mast  and 
a  single  antenna,  is  of  a  small  size  compared  with  that  required  to  represent  an  actual 
mast'antenna  structure  However,  the  model  has  all  the  features  necessary  to  allow  the 
assessment  and  critical  analysis  of  the  component  mode  synthesis  methods 

Specifically  ,  the  three  synthesis  methods  will  each  be  used  in  the  following  analy  ses 

(1)  Calculation  of  mast/antenna  coupled  system  modal  parameters:  This  is  the 
fundamental  assessment  of  a  method's  accuracy.  Prior  to  performing  the  synthesis,  modal 
parameters  are  calculated  for  the  antenna  model  and  the  mast  model  The  appropriate 
component  representation  is  generated  and  the  mast/antenna  system  is  synthesized  The 
coupled  system  natural  frequencies  are  calculated  and  are  compared  with  the  natural 
frequencies  calculated  using  a  standard  finite  element  procedure  The  standard  finite 
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element  procedure  means  the  assembly  of  a  single  model  representing  the  total 
mast/antenna  system  A  comparison  of  floating  point  operations  (FLOPS)  accumulated  in 
all  cases  is  also  provided  This  comparison  will  demonstrate  the  computational  advantage 
of  the  synthesis  methods,  an  advantage  critical  to  the  development  of  an  automated 
design  analysis  system 

Using  the  synthesized  mast 'antenna  model,  the  following  analyses  are  presented 

(2)  Calculation  of  antenna  peak  displacement  due  to  harmonic  forcing:  A  simple 
harmonic  forcing  function  is  applied  to  the  mast  and  the  peak  displacement  of  the  antenna 
free  end  (’’tip”)  is  calculated,  again  using  all  three  component  mode  synthesis  methods,  as 
well  as  using  a  standard  finite  element  procedure 

(3)  Calculation  of  mast/antenna  interface  internal  stresses  due  to  harmonic 
forcing:  A  simple  harmonic  forcing  function  is  applied  to  the  mast  and  the  bending 
moment  and  shear  loads  in  the  mast'antenna  connection  are  calculated  Note  that  these 
internal  loads  are  directly  proportional  to  stress,  and  hence  are  the  critical  quantities  which 
must  be  calculated  in  order  to  assess  structural  survivability  These  calculations  are 
repeated  for  all  three  synthesis  formulations,  as  well  as  for  the  standard  finite  element 
procedure 
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II.  FORMULATION  OF  FINITE  ELEMENT  MODEL  AND  GENERAL 
COMPONENT  COUPLING  PROCEDURES 


As  discussed  in  the  Introduction,  the  finite  element  (FE)  procedure  will  be  employed 
to  generate  mathematical  models  of  the  components  (substructures)  involved,  namely  the 
mast  and  the  antenna  The  FE  procedure  produces  stiffness,  mass,  and  less  commonly, 
damping  matrices  which  represent  the  structural  dynamics  of  each  component  In  order  to 
faithfully  capture  the  geometric  and  material  complexities  of  these  components,  the  finite 
element  discretization  must  necessarily  involve  many  degrees-of-freedom  (DOF),  and 
hence  the  above  mentioned  system  matrices  can  be  quite  large  The  time  and  cost 
associated  with  the  extraction  of  the  modal  parameters  (natural  frequencies,  mode  shapes, 
and  modal  mass)  from  these  large  matrices  precludes  the  performance  of  the  repeated 
design  analyses  required  to  arrive  at  an  optimal  design  The  component  mode  synthesis 
methods  bypass  the  repeated  extraction  of  the  modal  parameters  for  a  complete 
mast 'antennae  system  by  directly  using  the  modal  parameters  calculated  for  each 
component  The  calculation  for  the  component  “modes"  is  performed  once  for  each 
component,  and  the  total  system  dynamics  are  synthesized  using  the  various  sets  of  modes 
so  calculated  The  synthesis  methods  not  only  provide  very  accurate  predictions  of 
dynamic  response,  but  also  provide  a  substantial  decrease  in  the  time  required  to  compute 
dynamic  response,  hence  allowing  the  performance  of  additional  design  analysis  iterations 


A.  FINITE  ELEMENT  FORMULATION 

Although  FE  modeling  typically  involves  the  full  range  of  element  types  available  (e  g 
beam,  plate,  shell),  for  purposes  of  this  report  the  antenna  and  the  mast  will  each  be 
modeled  using  beam  elements  only  This  model,  although  simple,  is  all  that  is  necessary  to 
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investigate  the  various  component  coupling  procedures  All  methods  presented  herein  are 
applicable  to  any  structural  model,  and  the  results  and  conclusions  presented  are  directly 
applicable  to  the  analysis  of  structural  systems  of  any  complexity 

Traditionally,  the  mast  and  antennae  are  modeled  together  as  a  system  Alternatively, 
the  mast  and  antennae  can  be  modeled  separately  By  modeling  the  mast  and  antennae 
separately,  several  benefits  arise 

•  Masts  and  antennae  are  generally  fabricated  by  different  defense  contractors 
Therefore,  modeling  the  mast  and  antennae  separately  would  best  pres?  this 
independence 

•  Modeling  the  mast  and  antennae  separately  would  permit  the  development  of  a 
single  data  file  containing  only  mast  design  specifications,  and  several  separate  data 
files  containing  antenna  design  specifications,  one  datafile  for  each  antenna  With 
this  modular,  component-based  approach  comes  the  flexibility  of  exchangin, 
antennae  and  or  changing  antenna  placement  This  allows  the  rapid  assessment  of 
many  mast  and  antenna  configurations  for  dynamic  response  characteristics 

•  By  modeling  the  mast  and  antennae  separately,  the  computational  efficiency 
increases  as  compared  to  modeling  the  mast  and  antennae  together  This 
computational  advantage  is  due  to  the  fact  that  the  cost  associated  with  the 
calculation  of  the  modal  parameters  for  a  single  structural  model  is  proportional  to 
the  cube  of  the  number  of  DOF  of  the  model  [3]  The  calculations  performed 
herein  demonstrate  this  comparison  between  a  total  mast/antenna  model  and  a 
model  derived  from  the  synthesis  of  mast  and  antenna  substructure  models  The 
benefit  is  associated  not  just  with  the  calculation  of  the  modal  parameters,  but  also 
with  the  calculation  of  dynamic  response  to  assess  UNDEX  survivability. 
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B.  GENERAL  COMPONENT  COUPLING  PROCEDURES 


The  term  “component  mode  synthesis"  refers  to  the  manner  in  which  each  substructure 
is  mathematically  represented  prior  to  coupling,  and  is  based  on  a  truncated  modal 
expansion  This  representation  is  most  familiar  in  the  context  of  the  calculation  of  dynamic 
response  Here,  the  dynamic  response  of  a  structure  can  be  written  as  a  linear  combination 
of  the  mode  shapes  calculated  for  the  structure  If  the  frequency  range  of  excitation  is 
contained  in  the  frequency  range  of  the  calculated  modes,  then  the  dynamic  response 
calculated  using  the  modes  will  be  of  acceptable  accuracy  Of  course,  the  question  of  how 
man)  modes  to  retain  is  non-trivial  and  problem  specific  However,  the  computational 
efficiency  of  a  modal  approach  to  structural  dynamics  including  the  component  mode 
synthesis  methods  to  be  presented,  comes  from  the  retention  of  a  number  of  modes  w  hich 
constitute  a  mathematical  model  much  smaller  than  the  original  mass  and  stiffness  matrices 
from  which  the  modes  were  calculated 

Component  mode  synthesis  makes  use  of  several  types  of  vibrational  mode  shapes, 
distinguished  by  the  boundary  conditions  imposed  on  the  substructure  prior  to  the 
calculation  of  these  mode  shapes  In  addition  to  these  vibration  mode  shapes,  the  various 
component  mode  synthesis  methods  require  additional  types  of  mode  shapes  to  be 
calculated  and  included  with  the  vibration  modes  Therefore,  the  term  "component  mode 
synthesis"  (CMS)  is  a  suitable  name  a  single  structure  is  synthesized  from  separate 
substructures  and  each  substructure  is  mathematically  represented  by  an  appropriate  set  of 
mode  shapes  calculated  from  the  finite  element  model  of  each  substructure.  The  follow  ing 
are  definitions  of  the  various  types  of  mode  shapes  that  are  used  in  the  component  mode 
synthesis  formulations  investigated  herein 
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1.  FREE  INTERFACE  NORMAL  MODES 


The  free  interface  norma!  modes  are  the  modes  of  the  component  when 
unrestrained  at  all  interface  degrees  of  freedom  The  interface  coordinates  are  denoted  b> 
the  subscript  "I"  and  the  interna!  coordinates  are  denoted  by  the  subscript  "O"  The 
interface  coordinates  "I"  are  the  coordinates  where  the  substructures  are  coupled  The 
internal  coordinates  are  all  coordinates  that  are  not  interface  coordinates  Free  interface 
normal  modes  are  calculated  by  solving  the  following  eigenvalue  problem 

fK-XM]{<l>N}={0}  (1) 

The  stiffness  and  mass  matrices  in  Eq  (1 )  are  partitioned  as  follows 


The  number  of  equations  defined  by  Eq  (1 )  is  equal  to  the  number  of  row  s  or  columns 
in  the  mass  and  stiffness  matrices  The  number  of  columns  or  rows  in  [K]  or  [M]  equals 
the  number  of  DOF  of  the  component  in  physical  coordinates 

2.  FIXED  INTERFACE  NORMAL  MODES 

The  fixed  interface  modes  are  the  modes  of  the  component  restrained  at  its 
interface  degrees  of  freedom  The  fixed  interface  normal  modes  have  the  follow  ing  form 
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The  mode  shapes  in  the  upper  partition  of  the  matrix  of  fixed  interface  normal  modes,  or 
{*'}.  are  obtained  from  the  solution  to  the  following  eigenvalue  problem 

[K00-XM00  ]{*'}={()}  (2) 

In  words,  the  matrix  of  fixed  interface  normal  modes  is  a  partitioned  matrix  consisting 
of  the  mode  shapes  obtained  from  the  solution  of  Eq  (2)  in  the  upper  partition,  and  a 
matrix  of  zeros  in  the  lower  partition  The  zeros  mean  zero  displacement  at  the  interface 
The  number  of  rows  of  the  matrix  of  zeros  is  equal  to  the  number  of  interface  coordinates, 
while  the  number  of  columns  is  equal  to  the  number  of  internal  coordinates  of  the 
substructure 

Both  the  free  interface  norma!  modes  and  the  modes  obtained  in  the  solution  to  Eq  (2) 
are  unity  modal  mass  normalized  such  that  the  following  property  is  satisfied 


(3) 

(4) 

3.  STATIC  CONSTRAINT  MODES 

Static  constraint  modes  are  calculated  by  enforcing  a  unit  deflection  on  each 
interface  DOF  while  holding  all  other  DOF  restrained  Calculating  the  resulting 
displacements  of  the  interna!  coordinates  defines  the  static  constraint  modes  If  it  can  be 
assumed  that  no  external  forces  or  inertial  forces  are  applied  to  the  internal  degrees  of 
freedom,  as  in  a  static's  problem,  the  shapes,  or  [<!>'  ],  can  be  calculated  from  the  stiffness 
matrix  as  follows 


-K 


-i  K 

oo^oi  (  > 

I  F'* 


(5) 
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4.  RIGID  BODY  MODES 

Rigid  body  modes  are  associated  with  systems  that  are  not  constrained  Rigid 
body  modes  have  zero  frequency  They  can  be  solved  for  using  equation  (2 )  They  can 
also  be  solved  for  in  a  manner  identical  as  that  for  static  constraint  modes,  provided  that 
the  number  of  coordinates  retained  is  equal  to  the  number  of  rigid  body  modes  For 
purposes  this  report,  rigid  body  modes  will  be  solv  ed  for  in  accordance  with  equation  (5) 

5.  RESIDl'AL  FLEXIBILITY  MODES 

Before  defining  "residua!"  flexibility,  the  concept  of  flexibility  must  first  be  defined 
The  flexibility  of  a  restrained  structure  (i  e  a  structure  whose  stiffness  matrix  is  of  full 
rank)  is  the  inverse  of  the  stiffness  of  the  structure  By  inverting  the  stiffness  matrix,  one 
obtains  the  flexibility  matrix  as  follows 


[G]=[K]';  (6) 

Equation  (6)  can  also  be  written  as  follows 

(G]=[4>'I  A  ]■’[*' f  (7 > 

The  residual  flexibility  matrix  is  obtained  from  the  flexibility  matrix,  the  kept  free  interface 
normal  modes,  and  the  inverse  of  the  natural  frequencies  as  follow  s 

[Grf]  =[G] —[$,'][  Ak  ]"[4>ks]1  =[*„'][  K  ]"[*osr  (8) 

The  residual  flexibility  modes  are  the  portion  of  the  exact  static  flexibility  shapes  that  are 
not  represented  by  a  set  of  retained  modes  Residual  flexibility  modes  require  the 
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knowledge  of  other  modes  that  are  retained  in  the  model  and  are  dependent  upon  the 
retained  modes  There  are  two  ways  to  calculate  the  residual  flexibility  modes 

1)  If  the  structure  is  grounded,  such  as  the  mast,  then  the  stiffness  matrix  is  full  rank 


and  invertible 


By  post-multiplying  Eq  (8)  by 


one  obtains  the  residual  flexibility 


modes  for  restrained  substructures 


W  =[Grf  j 

L1  u  . 


(9) 


2)  However,  If  the  component  is  not  grounded  before  assembly,  such  as  an  antenna, 
then  an  inertia  relief  solution  must  be  calculated  to  determine  the  flexibility  matrix  as 
follows 


[G]  =[I-*r*rm]T[K*][i-¥r*k’m]  (10) 

w  here  are  the  rigid  body  modes  of  the  structure 

[K  *]  is  formed  by  inverting  the  restrained  or  internal  degrees  of  freedom  in  the  system  in 
the  following  way 


[K  *] 


0  0 

?  *C. 


(11) 


This  "new  "  [G]  or  flexibility  matrix  in  Eq  (10)  is  free  of  rigid  body  modes  The  Craig- 
Chang  formulation,  which  will  be  presented  in  the  next  section  uses  free-interface  normal 
modes  and  residual  flexibility  modes  The  residual  flexibility  modes  of  unrestrained 
substructures  are  obtained  from  the  neglected  or  deleted  free  interface  normal  modes  just 
like  they  were  obtained  from  a  substructure  that  is  restrained  The  only  difference  is  that 
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the  flexibility  matrix  obtained  in  Eq  (1 0)  is  used  Residual  flexibility  modes  are  calculated 
by  computing  the  static  flexibility  and  subtracting  the  flexibility  due  to  the  retained  modes 
The  residual  flexibility  modes  are  obtained  from  the  flexibility  matrix  in  Eq  (10)  in  the 
same  way  that  they  were  obtained  in  Eq  (9) 


W  =[Grf] 


(0, 


■n  j 


(12) 
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III.  COMPONENT  MODE  SYNTHESIS  FORMI  CATIONS 


A.  CRAIG-BAMPTON  FORM  ELATION 

There  are  three  substructure  coupling  procedures  that  serve  as  potential  candidates  to 
be  used  in  the  mast 'antenna  synthesis  The  Craig-Bampton  formulation,  the  Craig-Chang 
residual  flexibility  method,  and  the  MacNeal  residual  flexibility  method  In  this  section, 
the  Craig-Bampton  formulation  will  be  presented,  while  in  Section  B  the  Craig-Chang  and 
MacNeal  residual  flexibility  methods  will  be  presented  together  because  of  the  similarities 
in  the  methods 

The  Craig-Bampton  reduction  procedure  uses  a  combination  of  static  constraint 
modes  and  fixed  interface  normal  modes  to  represent  the  component  model  Both  the 
static  constraint  modes  and  the  fixed  interface  normal  modes  are  obtained  from  the  finite 
element  substructure  models  This  combined  set  of  mode  shapes  will  be  used  to  transform 
the  original  large  order  substructure  mass  and  stiffness  matrices  down  to  a  significantly 
smaller  size,  a  size  equal  to  the  number  of  mode  shapes  included  in  the  transformation 
matrix  The  transformation  matrix.  M  for  the  Craig-Bampton  formulation  contains  the 
shape  functions  as  its  columns  as  follows 


This  transformation  matrix  is  obtained  for  each  substructure  in  the  system  The  size  of 
the  static  constraint  mode  partition  of  the  transformation  matrix  is  always  fixed,  because 
the  number  of  columns  corresponds  to  the  number  of  interface  degrees  of  freedom 
However,  the  fixed  interface  normal  mode  partition  is  not  held  constant  The  number  of 
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columns  can  range  as  low  as  one  column  if  only  one  fixed  interface  norma!  mode  is 
retained,  or  as  high  "m"  columns  where  "m"  is  the  total  number  of  internal  degrees  of 
freedom  The  size  of  the  transformation  matrix  depends  upon  how  many  modes  are 
required  to  accurately  represent  the  physical  dynamic  response  of  the  system  when 
subjected  to  a  forced  input  Retaining  fewer  modes  than  the  total  possible  modes  available 
is  referred  to  as  "modal  truncation,"  and  provides  the  computational  efficiency  of  the 
method  Retaining  fewer  modes  than  the  total  amount  of  modes  available  means  fewer 
calculations  required  in  conducting  the  dynamic  analysis  On  the  other  hand,  if  the  number 
of  modes  retained  are  not  sufficient  to  accurately  determine  the  dynamic  response,  then 
the  benefits  of  reduced  compute  times  do  not  outw  eigh  the  magnitude  of  error  obtained  in 
the  analysis  Therefore,  while  the  benefits  of  modal  truncation  are  important  in  shortening 
compute  times,  they  are  not  as  important  as  obtaining  accurate  results  In  terms  of 
computational  efficiency,  large  benefits  can  be  achieved  using  this  method  if  only  the 
lower  range  of  frequencies  is  of  interest  This  is  applicable  to  the  mast  which  is  subjected 
to  typically  low  forcing  frequencies  By  retaining  a  few  of  each  of  the  component  modes, 
an  accurate  assessment  of  the  dynamic  response  of  the  mast  and  antenna  is  obtained  The 
examples  contained  in  section  IN'  and  V  demonstrate  how  by  retaining  just  a  few  modes  of 
each  substructure,  accurate  results  can  be  obtained 

By  pre-  and  post-multiplying  the  respective  substructure  mass  and  stiffness  matrices  by 
the  transformation  matrix  in  the  following  manner,  one  obtains  the  reduced  component 
model  as  follows 


(14a) 

(14b) 
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The  partitions  of  the  reduced  mass  and  stiffness  matrices  are  expressed  as  follows 


=[A  J 

k;  =k"  =[o] 
k;  =[ku  -k1ukc„,-kui] 

m;  =Dkk] 

M;  =M;’ 

The  term  reduced,  designated  by  the  subscript  V,  means  that  the  resulting  mass  and 
stiffness  matrices.  Ik.]  &  [m,]  are  of  smaller  dimension  than  the  original  matrix 
Although  this  transformation  matrix  reduces  the  size  of  the  component  model,  it  does  not 
assemble  the  individual  substructure  models  There  is  a  second  transformation  matrix  that 
synthesizes  the  substructures  to  produce  the  total  system  by  enforcing  the  compatibility 
and  equilibrium  of  the  interface  coordinates  as  follows 

fi'Mv)  05a) 

{r;}--{F;}  <i5w 

The  compatibility  of  interface  coordinates  denoted  by  £q(15a)  implies  that  the 
displacement  at  the  interface  of  structure  1  equals  the  displacement  at  the  interface  of 
structure  2  Likewise,  the  equilibrium  of  interface  forces  denoted  by  Eq  (15b)  implies  that 
the  sum  of  the  forces  at  the  interface  are  equal  but  in  acting  in  opposite  directions  Since 
the  static  constraint  modes  are  independent  of  the  fixed  interface  normal  modes,  the 
transformation  matrix,  or  [T,  ].  takes  the  following  form: 
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where  superscript  "1"  refers  to  substructure  1  (i  e ,  mast),  and  superscript  "2”  refers  to 
substructure  2  (i  e  ,  antenna)  In  Eq  (16).  p0  represents  the  physical  internal  coordinates, 
and  p,  represents  the  physical  interface  coordinates 

The  system  mode!  is  obtained  from  an  uncoupled  mass  and  stiffness  matrix  These 
uncoupled  mass  and  stiffness  matrices  are  themselves  formed  from  the  reduced  mass  and 
stiffness  matrices  front  each  substructure  as  follows: 


[K.]  = 


k;  o 

0  K; 


(17a,b) 


where  the  subscript  “u"  denotes  “uncoupled  “ 

The  coupled  system  mass  and  stiffness  matrices  are  obtained  by  pre-multiplying  the 
uncoupled  mass  and  stiffness  matrices  in  Eq  (17a,b)  by  the  transpose  of  the 
transformation  matrix  in  Eq  (16)  and  then  post-multiplying  the  uncoupled  mass  and 
stiffness  matrices  b\  the  transformation  matrix. 


(18a.b) 


where  the  subscript  “s"  denotes  “system.” 
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The  resulting  coupled  system  mass  and  stiffness  matrices  are  of  the  form 
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(19a,b) 


where  "1"  represents  the  interface  coordinates,  and  "K"  represents  the  kept  internal 
coordinates  [4]  The  respective  partitions  of  the  system  mass  and  stiffness  matrices  are 
expressed  as  follows 

ML  =[lL] 

ml  =DL] 

ML,  =m»  =[+','(ML(-K;:„K„t)’  +ml,i] 

ML,  -ML’  =hV'(ML(-Ki,K0,)-’  +M;„)] 

m„  -kL.Kd,  -k;L.k„,)’  +m;„)+mU-KooK0,V  +m;,  + 

(-Kn'nK0,),’(ML,<-KiK„l)!  +ML,)+ML(-KL,K0,)!  +Mj] 

kl  =[a  :,,] 

KL  =[a  L] 

=[K|,  -kl,kc:„'k;„  +KL  -KL.KLLKL,] 

Theie  are  a  number  of  advantages  to  the  Craig-Bampton  component  mode 
representation  The  first,  which  is  especially  beneficial  to  the  analysis  of  the  mast/antenna 
system,  is  that  the  reduced  DOF  system  contain  the  interface  DOF  explicitly.  This  makes 
it  very  easy  to  couple  mast  and  antenna  substructures  In  the  figure  on  the  following  page 
is  an  illustration  of  the  mast/antenna  system  used  for  the  examples  in  this  report  Along 
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the  length  of  the  cross  bar  are  various  node  positions  These  node  positions  serve  to 
connect  the  beams  that  represent  the  cross  bar,  and  can  also  serve  as  nodes  to  connect 
various  antennae  to  the  crossbar  By  specify  ing  different  "connection  coordinates”  (i  e , 
the  "I"  coordinates),  and  with  separate  mast  and  antenna  data  files  contained  in  the  Jibraty, 
the  engineer  can  quickly  couple  various  antennae  with  the  mast  and  rapidly  determine 
dynamic  response  Should  the  location  of  the  antenna  placement  not  be  suitable,  the 
engineer  can  specify  a  new  set  of  interface  coordinates  along  the  cross  bar,  plug  in  the 
antenna  at  the  new  location,  and  rapidly  calculate  a  system  from  which  a  new  dynamic 
response  can  be  calculated 

Because  of  these  advantages,  and  because  the  Craig-Bampton  component  mode 
representation  tends  to  result  in  accurate  system  frequencies,  as  will  be  shown,  this  is  a 
widely  used  method  Additionally,  the  NASTRAX  superelement  scheme  uses  the  Craig- 
Bampton  component  mode  representation  with  minor  extensions  as  a  solution  path  to  the 
dynamic  response  problem 

B.  CRA1G-CHANG  AND  MACNEAL  RESIDUAL  FLEXIBILITY 

FORMULATIONS 

The  Craig-Chang  and  MacNeal  residua!  flexibility  formulations  will  now  be  discussed 
Due  to  the  similarity  in  the  methods,  the  Craig-Chang  procedure  will  be  presented  first, 
and  the  modification  of  this  method  to  produce  the  final  system  of  equations  of  the 
MacNeal  method  will  be  discussed  subsequently  While  the  Craig-Bampton  representation 
uses  a  combination  of  static  constraint  modes  and  fixed  interface  normal  modes,  the  Craig- 
Chang  residual  flexibility  formulation  combines  free  interface  mode  shapes  with  residual 
flexibility  shapes,  thus  the  name  residua!  flexibility  method  The  transformation  matrix 
which  is  used  to  reduce  the  component  mass  and  stiffness  matrices  contains  columns  of 
the  retained  or  kept  free  interface  normal  modes  and  residual  flexibility  modes 
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The  transformation  matrix,  or  [T,l  is  shown  as  follow  s 


(20) 


In  the  same  way  that  we  reduced  the  Craig-Bampton  components,  we  reduce  the 
Craig-Chang  components  except  now  using  the  transformation  matrix  of  Eq  (20)  Again, 
each  component  has  its  own  transformation  matrix  A  reduction  in  component  matrix  size 
is  achieved  by  retaining  few  er  than  the  total  number  of  free  interface  normal  modes  The 
number  of  residual  flexibility  mode  shapes  is  fixed,  and  equals  the  number  of  interface 
degrees  of  freedom 

Just  like  the  Craig-Bampton  formulation,  the  purpose  of  the  transformation  matrix  is 
to  reduced  the  respective  component  models  before  synthesis  of  the  system  model  This 
transformation  matrix  does  not  synthesize  the  substructures  Another  transformation 
matrix  is  employed  to  synthesize  the  substructures  As  in  the  Craig-Bampton  formulation, 
this  second  transformation  matrix  results  from  satisfying  compatibility  and  equilibrium 
equations 


(21a) 

{f;  >=  -{F,! }  (21b) 


I  niike  the  Craig-Bampton  component  mode  representation  where  the  static  constraint 
modes  .ire  independent  of  the  fixed  interface  normal  modes,  the  residua!  flexibility  modes 
are  dependent  upon  the  free  interface  normal  modes,  and  a  simple  boolean  matrix  will  not 
synthesize  the  substructures  Although  not  derived  in  this  report,  the  second 
transformation  matrix,  or  [T,].  is  as  follow  s 
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[T]  = 


-k  4>' 

~  l*kl 

k  d>2 
*1*KI 

-ML 

HI 

[o] 

[0] 

[i] 

(22) 


where  k,  =('i,1'1 

As  in  the  Craig-Bampton  procedure,  the  uncoupled  mass  and  stiffness  matrices  are 
assembled  from  the  reduced  component  mass  and  stiffness  matrices 


M»D 
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(23a,b) 


By  pre-  and  post-multiplying  both  the  uncoupled  mass  matrix  represented  by  Eq  (23a) 
and  the  uncoupled  stiffness  matrix  represented  by  Eq  (23b)  by  [T,]t  and  [T,] 
respectively,  in  the  same  manner  that  the  Craig-Bampton  component  mode  representation 
w  as  coupled,  the  system  equations  of  motion  are  obtained  as  follow  s 


The  partitions  of  the  system  mass  and  stiffness  matrices  represented  by  Eq  (24a, b)  are 
given  as  follows: 

M„  =[4  +*£n>,*L]  K„  =[A  ^ 

M„  - Mj,  =H>£m,iL]  K„  -Kj,  =[-<k,^] 

=[lL  +*kfn,H]  K„  =[A^  4+^,*’*] 
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M„  -  Mil  =H>£"',*i;] 

U.  m  1^  u,  ] 


Kn  =[A  L  +*^i*1k] 
K,i  =Kl,  =[-4>ik,*^] 

K„  =[AL  +<k,«4] 


where  m(  =k,(Mj)n  +M^n)k,  [5]  The  inertia  due  to  high  order  free  interface  normal 
modes  is  represented  by  “nv  ” 

As  stated  in  the  beginning  of  this  section,  the  MacNeal  component  mode 
representation  would  be  presented  By  neglecting  the  inertia  due  to  high-order  free- 
interface  normal  modes,  (i  e  “m  "),  one  obtains  the  MacNeal  mass  and  stiffness  matrices 
as  follow  s 

M„  =[l’J  Kn  “[A  L 

M[:  =Mj,  =[0]  K);  =K],  =[-4>£k,«^] 

M:.’  “Dkx]  Kj2  =[A  k  +*ikk1*L] 

As  will  be  demonstrated  through  the  examples,  the  effect  of  neglecting  “m/‘  is 
important  when  predicting  the  higher  frequencies  [6]  The  MacNeal  representation  is 
accurate  in  the  lower  and  mid  frequency  range,  but  less  accurate  in  the  higher  frequency 
range 

Just  like  the  C'raig-Bampton  component  mode  representations  there  are  several 
advantages  to  both  residual  flexibility  methods  By  analyzing  both  the  Craig-Chang  and 
MacNeal  system  mass  and  stiffness  matrices,  one  can  see  that  the  final  system  coordinates 
are  just  the  free  interface  normal  mode  coordinates  from  each  substructure  This  resulted 
from  the  operations  that  were  conducted  in  forming  the  system  mass  and  stiffness 
matrices.  Additionally,  since  the  residual  flexibility  modes  account  for  the  static  flexibility 
of  all  modes,  the  methods  are  statically  exact  The  procedure  is  applicable  to  the  mast 
and  antenna  problem  Since  the  connection  coordinates  are  not  explicitly  retained,  this 
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connection  points,  and  the  residual  flexibility  method  does  not  retain  the  connection 
coordinates  explicitly,  a  reduction  in  compute  times  result,  an  advantage  not  found  in 
other  substructure  coupling  procedures 


26 


IV.  NUMERICAL  VERIFICATION 


In  this  section  some  numerical  examples  which  demonstrate  the  use  of  the  three 
component  mode  synthesis  procedures  will  be  presented  A  "mock-up”  mast  and  antenna 
system  consisting  of  1 7  elements  was  assembled  using  a  standard  finite  element  procedure 
and  all  three  component  mode  synthesis  procedures  Although  the  mast/antenna  system 
that  is  used  in  the  examples  is  modeled  with  just  a  few  elements,  the  resulting  models  are 
large  enough  to  allow  the  effects  of  mode  truncation  to  be  assessed  Again,  it  is  not  the 
intent  of  this  report  to  solve  a  base  excitation  problem  on  a  realistic  mast  and  antenna 
model,  but  rather  to  demonstrate  how  CMS  can  be  used  when  performing  dynamic 
analyses  for  design  purposes 

Example  1  compares  natural  frequency  calculations  for  the  total  mast.'antenna  system 
as  computed  using  the  three  CMS  formulations  as  well  as  using  a  standard  FE  procedure 
Tables  1-3  contain  the  results  of  this  comparison  In  the  tables,  each  row  contains  the 
estimate  of  a  mode  frequency  The  first  column  contains  the  mast/antenna  system  natural 
frequency  estimates  as  calculated  using  the  standard  FE  procedure,  and  serves  as  the 
reference  value  against  which  the  CMS  natural  frequency  estimates  are  to  be  compared 
Columns  2  through  4  contain  the  analogous  natural  frequency  estimates  and  percent  error, 
as  calculated  from  the  mast/antenna  system  model  synthesized  using  each  of  the  three 
CMS  procedures  Also  included  in  the  column  headings  are  floating  point  operations 
(FLOPS)  counts  which  provide  a  measure  of  the  number  of  calculations  required  to 
assemble  the  mast/antenna  system  and  calculate  the  natural  frequencies  and  mode  shapes 

Table  1  presents  the  system  frequency  comparison  where  18%  of  the  available  mast 
modes  are  retained,  and  22%  of  the  available  antenna  modes  are  retained.  Table  2  repeats 
the  calculations  with  42%  of  the  available  mast  modes  retained  and  39%  of  the  available 
antenna  modes  retained,  and  Table  3  repeats  the  calculations  with  79%  of  the  available 
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mast  modes  retained  and  67%  of  the  available  antenna  modes  retained  Note  that  each 
subsequent  table  presents  comparisons  for  an  increasing  number  of  mode  frequencies  due 
to  the  fact  that  an  increase  in  the  number  of  retained  component  modes  makes  possible  an 
increase  in  the  number  of  system  modes  which  may  be  calculated 

Note  that  in  the  FEM  model,  the  FLOPS  count  stays  fixed  at  slightly  over  2-10fc 
This  is  a  rather  small  number  as  the  model  is  a  small  model  when  compared  to  one  that  a 
design  engineer  would  generate  for  analysis  of  an  actual  mast/antenna  assembly  It  is 
noteworthy  to  state  that  in  this  particular  model.  07  -106  FLOPS  were  expended  in 
computing  the  combination  of  fixed  interface  normal  modes  and  static  constraint  modes 
using  the  Craig-Bampton  procedure  Additionally,  0  8-10*  FLOPS  were  expended  in 
computing  the  free  interface  normal  modes  using  the  Craig-Chang  and  MacNeal 
procedures  Theoretically,  once  the  various  vibrational  modes  have  been  found,  they  need 
not  be  calculated  again  Note  also  that  in  all  three  models,  these  figures  comprise  a 
significant  portion  of  the  total  FLOPS 

From  Tables  1  through  3,  it  is  seen  that  all  three  methods  produce  excellent  frequency 
predictions  All  three  methods  demonstrate  sudden  increases  in  frequency  error  above  a 
certain  mode  This  reveals  the  extent  to  which  the  retained  component  modes  accurately 
represent  the  dynamics  of  the  synthesized  mast/antenna  system  In  the  MacNeal 
procedure,  the  percentage  error  exceeded  100%  when  calculating  the  highest  mode  This 
error  in  predicting  the  highest  frequency  mode  could  possibly  be  attributed  to  neglecting 
the  inertia  due  to  high-order  free  interface  normal  modes  Note  that  by  neglecting  the 
inertia,  that  the  accuracy  in  predicting  natural  frequencies  is  only  effected  at  the  last  few 
modes 

As  stated  in  the  previous  paragraph,  all  three  methods  produced  excellent  results  in 
predicting  natural  frequencies,  but  with  less  cost  in  terms  of  number  of  computations  as 
compared  to  the  standard  FE  calculation.  The  Craig-Chang  procedure  in  general  provided 
the  greatest  number  of  natural  frequencies  with  error  less  than  or  equal  to  0  1%  (in  Hz) 
(see  figure  2)  However,  the  Craig-Bampton  procedure  yielded  the  same  number  of 
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frequencies  with  error  less  than  or  equal  to  0  1%  as  the  Craig-Chang  procedure  when 
retaining  a  large  of  number  of  component  modes,  but  at  a  slightly  more  cost  than  the 
Craig-Chang  procedure 


Table  1.  Coupled  System  Natural  Frequency  Comparison 
Percent  of  Available  Mast  Modes  Retained:  18% 
Percent  of  Available  Antenna  Modes  Retained:  22% 
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Table  3.  Coupled  System  Natural  Frequency  Comparison 
Percent  of  Available  Mast  Modes  Retained:  79% 
Percent  of  Available  Antenna  Modes  Retained:  67% 
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V.  BASE  EXCITATION  FORMULATIONS 


In  the  previous  section,  the  natural  frequencies  of  the  mast  and  antenna  system  were 
calculated  for  increasing  number  of  retained  component  modes  Natural  frequencies  and 
mode  shapes  are  important  modal  parameters  and  are  fundamental  in  solving  for  the 
forced  response  of  a  system  As  demonstrated  in  the  previous  section,  accurate  natural 
frequencies  of  a  system  can  be  obtained  using  CMS  at  a  cost  less  than  that  associated  with 
standard  FE  modeling 

In  this  section,  two  base  excitation  formulations  will  be  presented  The  first 
formulation  requires  the  knowledge  of  the  acceleration  of  the  mast  base  coordinates  (i  e 
the  coordinates  where  the  mast  and  ship  are  coupled)  as  a  function  of  time  In  other 
words,  the  formulation  requires  that  the  acceleration  time  history  of  the  base  coordinates 
be  know  n  The  second  formulation  requires  the  knowledge  of  the  displacement  of  the 
base  coordinates  as  a  function  of  time,  or  the  displacement  time  history'  of  the  base 
coordinates  Using  both  formulations,  numerical  convergence  assessments  will  be  made, 
and  the  benefits  that  CMS  has  to  offer  the  mast/antenna  design  process  will  be 
demonstrated 

A.  BASE  EXCITATION  FROM  PRESCRIBED  ACCELERATION 

Once  the  FE  program  has  numerically  assembled  the  mast  mass  and  stiffness  matrices, 
and  the  acceleration  of  the  base  coordinates  are  specified  as  a  function  of  time,  the  base 
excitation  problem  can  be  derived  from  the  following  equation  of  motion: 

[Moo  Mob 
IMbo  Muj, 


K«o  *t>  | 1 

Nllu  Kbb.  Xb  1  ®  J 


(25) 
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where  subscript  "O"  represents  the  interior  coordinates,  and  subscript  ”B"  represents  the 
base  coordinates 

Solving  the  top  row  of  equations  in  Eq  (25)  the  following  is  obtained 

MocX0  "F^OBXB  +K00X0  "^^OBXB  —  f()  (26) 

Since  the  acceleration  of  the  base  is  prescribed,  the  base  acceleration  term  will  be  moved 
to  the  right  hand  of  the  equals  sign  to  obtain  the  following 

Muu*o  +Kuoxo  "^^odxe  Fu  Mubxb  (27) 

From  the  bottom  row  ,  the  following  equation  is  obtained 

Mbuxo  +Mbl.xl.  +K-Bl>x1.  +KbBxb  =0  (28) 

From  Eq  (28),  the  following  relation  is  obtained  for  the  base  displacement 

xp,  —  [-^roxo  ■^^■nfx()]  (29) 

Equation  (29)  is  now  substituted  into  Eq  (27).  and  after  simplifying,  the  following 
equation  of  motion  in  terms  of  the  interior  coordinates  is  obtained  as  follow  s 

[Moo  “K0|,K.nr,M|)o  +[k()(,  "“Kyj.KjjpKjj,, Fc  +[kncKDnMur  (->0) 

Equation  (30)  is  the  system  equation  of  motion  of  the  internal  coordinates  in  terms  of  the 
prescribed  base  acceleration 
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The  physical  coordinates  are  transformed  into  modal  coordinates  by  the  following  relation 

{XoH+'l  {l«}  <3I> 

where  {xo]  represents  the  vector  of  physical  internal  coordinates,  {q0}  represents  the 
vector  of  modal  internal  coordinates,  and  [$s]  represents  the  matrix  of  unity  modal  mass 
free  interface  normal  modes 

Equation  (30)  is  then  pre-multiplied  by  the  transpose  of  the  matrix  of  normal  modes  to 
obtain  the  follow  ing  modal  system  of  equations 

il]{U+[ A  &„}-(/  !  (32) 

where  {qv, }  represents  the  vector  of  modal  accelerations  and  {/-'}  represents  the  vector  of 
modal  forces 

The  relation  in  Eq  (32)  is  solved  using  standard  modal  decomposition  techniques 
Once  the  modal  response  is  obtained.  Eq  (31)  is  used  to  obtain  the  response  in  physical 
coordinates 

In  what  follows,  the  relation  in  Eq  (30)  is  used  in  conjunction  with  a  standard  FE 
model  of  the  total  mast/antenna  system  to  perform  the  prescribed  base  acceleration 
dynamic  analysis  This  analysis  serv  es  as  the  reference  against  w  hich  the  results  of  various 
CMS  formulations  are  to  be  compared  In  each  CMS  formulation,  Eq.  (30)  is  used  to 
define  the  mast  component  model,  as  only  the  mast  has  prescribed  base  accelerations 
The  acceleration  at  the  base  is  a  prescribed  harmonic  input.  The  base  acceleration  is 
taken  as 


(33) 
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where  V"  represents  the  forcing  frequency  expressed  in  rad/sec,  and  {XH}  represents  the 
vector  of  amplitudes  of  the  base  displacement  expressed  in  inches 

In  the  following  section,  two  numerical  examples  are  provided  It  is  the  intent  of  the 
examples  to 

( 1 )  compare  the  results  obtained  from  the  three  CMS  formulations 

(2)  and  demonstrate  the  benefits  of  using  CMS  versus  standard  FE  modeling 
when  solving  base  excitation  problems 

1.  TIP  DEFLECTION  CALCULATION 

On  the  follow  ing  page  is  a  diagram  of  the  mast  and  antenna  system  being  subjected 
to  base  excitation  The  excitation  was  performed  at  two  different  and  arbitrarily  selected 
frequencies  The  first  frequency  was  at  8  95  Hz  This  frequency  falls  between  mode  I  and 
mode  2  of  the  total  mast/antenna  system  The  second  frequency  was  betw  een  mode  9  and 
mode  10  at  2 19  Hz  The  wide  spread  in  the  frequencies  was  intended  to  demonstrate  that 
many  more  modes  need  to  be  retained  when  calculating  the  response  to  higher  frequency 
excitation  as  compared  with  lower  frequency  excitations  In  this  example,  the  antenna  tip 
deflection  was  calculated  using  the  standard  FE  procedure  and  the  three  CMS  procedures 
The  percent  error  in  antenna  tip  deflection  was  plotted  versus  the  percent  of  available 
component  modes  retained  (see  figures  4-7)  The  calculations  were  performed  twice  In 
the  first  calculation,  mast  modes  were  truncated  while  retaining  all  of  the  available  antenna 
modes,  and  antenna  modes  were  truncated  in  the  second  calculation  while  retaining  all  of 
the  available  mast  modes  When  the  mast  was  subjected  to  the  forced  input  at  the  lower 
frequency  the  Craig-Bampton  and  Craig-Chang  procedure  >  aided  results  which 
converged  more  rapidly  to  the  exact  answer  than  the  MacNeal  procedure 
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Figure  3:  The  Coupled  Mast  and  Antenna  Subjected  to  Base 
Excitation  at  8.95  Hz 
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Although,  hard  to  determine  from  figures  4  and  5.  the  Craig-Bampton  procedure  yielded 
the  best  results  using  fewer  modes  than  the  Craig-Chang  procedure  in  both  the  mast 
truncation  and  antenna  truncation  runs  When  20%  of  the  available  mast  modes  were 
used,  all  three  procedures  predicted  a  tip  deflection  measurement  that  was  within  0  05% 
of  the  "exact"  value  (the  “exact"  value  was  calculated  using  the  standard  FE  procedure) 
This  accurate  assessment  was  obtained  at  a  cost  of  60%  of  the  number  of  calculations  that 
were  required  of  the  FE  solution  The  error  obtained  when  20%  of  the  available  mast 
modes  were  retained  was  significantly  small  Since  all  of  the  vibrational  mode  shapes  w  ere 
obtained  in  this  calculation  (i  e  fixed  interface  normal  modes,  free  interface  normal  modes, 
static  constraint  modes,  and  residual  flexibility  modes),  future  assessments  which  retain 
more  mast  modes  would  come  at  an  even  lesser  cost  than  the  initial  assessment 

When  the  mast  was  subjected  to  the  higher  forcing  frequency,  all  three  methods 
converged  more  slowly  as  compared  to  the  lower  forcing  frequency  when  truncating  both 
mast  modes  and  antenna  modes  (see  figures  6  and  7)  Since  the  forcing  frequency  was 
higher,  more  modes  needed  to  be  retained  in  order  to  obtain  accurate  results  From  the 
results  of  the  tip  deflection  calculations,  it  appears  that  the  combination  of  fixed  interface 
normal  modes  and  static  constraint  modes  have  led  to  the  higher  rate  of  convergence  using 
the  Craig-Bampton  procedure  However,  the  results  obtained  using  the  Craig-Chang 
procedure  compared  quite  well  with  the  results  obtained  using  the  Craig  Bampton  method 
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Figure  4:  Percent  error  in  antenna  tip  deflection  plotted  versus  the 

percent  of  available  mast  modes  retained.  (Forcing  frequency:  8.95  Hz) 
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Figure  5:  Percent  error  in  antenna  tip  deflection  plotted  versus  the 

percent  of  available  antenna  modes  retained.  (Forcing  Frequency:  8.95  Hz) 
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Figure  6:  Percent  error  in  antenna  tip  deflection  plotted  versus  the 
percent  of  available  mast  modes  retained.  (Forcing  Frequency:  219  Hz) 
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Figure  7:  Percent  error  in  antenna  tip  deflection  plotted  versus  the 
percent  of  available  antenna  modes  retained.  (Forcing  Frequency:  219  Hz) 
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2.  MOMENT  AND  SHEAR  CALCULATION 

In  order  to  assess  structural  survivability,  accurate  prediction  of  the  internal 
stresses  in  the  antenna  and  the  antenna'mast  interface  must  be  calculated  Therefore,  this 
section  demonstrates  the  calculation  of  the  internal  peak  dynamic  bending  moments  and 
shear  loads  in  the  antenna 

The  moment  and  shear  calculations  were  calculated  using  the  three  CMS  procedures 
and  the  results  are  compared  in  the  figures  The  percent  error  in  moment  and  shear  were 
plotted  versus  the  percent  of  available  mast  modes  retained  and  percent  of  available 
antenna  modes  retained  Again,  the  same  two  forcing  frequencies  used  in  the  tip 
deflection  calculation  are  used  here  in  the  moment  and  shear  calculations,  specifically  8  95 
and  219  Hz 

The  results  obtained  when  calculating  the  shear  and  moment  at  the  mast'antenna 
connection  mirror  the  results  of  the  tip  deflection  calculations  (see  figures  8-15)  Again, 
the  Craig-Bampton  procedure  yielded  results  that  converged  more  quickly  to  the  "exact" 
answer  (provided  by  standard  FE  calculations)  than  the  other  methods  However,  the 
results  obtained  using  the  Craig-Chang  procedure  were  quite  similar  to  those  obtained 
using  the  Craig-Bampton  method  Despite  a  large  initial  error  produced  by  the  MacNeal 
method  as  compared  to  the  other  two  methods,  nearly  "exact”  solutions  were  obtained  at 
a  cost  much  less  than  using  standard  FE  calculation  procedures  If  the  moment  and  shear 
at  the  mast'antenna  interface  exceeded  an  appropriate  failure  criteria,  the  antenna  can  be 
easily  relocated  from  the  end  node  to  another  node  along  the  cross  bar  by  redefining  the 
connection  coordinates  of  the  mast  New  moment  and  shear  calculations  would  be  made 
until  an  acceptable  response  obtained  Redefining  the  connection  coordinates,  synthesizing 
the  new  structure,  and  calculating  the  response  is  much  more  convenient  and 
computationally  efficient  than  reassembling  the  mast  and  antenna  system,  which  would  be 
required  using  standard  FE  procedures 


When  the  mast  and  antenna  system  were  subjected  to  the  forced  input  at  the  higher 
forcing  frequency,  the  rate  of  convergence  was  again  much  slower  than  that  which  was 
obtained  at  the  lower  excitation  frequency  However,  all  methods  yield  accurate  results 
at  a  computational  cost  less  than  using  the  standard  FE  procedure  with  the  higher  forcing 
frequency 


Percent  error  in  moment  at  mast  Percent  error  in  shear  at  matt  and 

anil  antenna  connection  antenna  connection 


Percent  error  in  shear  at  mast  and  antenna  connection  plotted  versus  the  percent  of 
available  mast  modes  retained 


1  60% 

4 

t  60% 

— a 

CRAJG-BAMPTON 

1  40% 

- 

-A  - 

craig-chang 

1  20% 

- 

1.00% 

- 

A 

MAONKAL 

0  80%, 

0  60% 

- 

0  40%, 

♦ 

A 

0  20% 

a  . _ 

0  00  Vo 

a 

~  -  -  -  f  -  ---  '  "g 

0% 

20% 

30%  40%  50% 

1 

o'* 

O 

<0 

70%  80% 

1 

Number  of  mast  modes  retained  (in  percent) 

Figure  8:  Percent  error  in  shear  at  mast  and  antenna  connection  plotted  versus 
the  percent  of  available  mast  modes  retained.  (Forcing  Frequency:  8.95  Hz) 


Percent  error  in  moment  at  mast  and  antenna  connection  plotted  sersus  the  percent 
of  as  ailable  mast  modes  retained 


12  00  V. 

* 

9  00 


6  00  -  - 


- •  —  CKA1G-BA\1PI<  '\ 

-  ♦  -  CRAlG-GHAW. 

a  M  ACM  At 


3.00*  - 


0.00% 


107c  20%  30%  40%  50%,  60%  70%  60% 

Number  of  mast  modes  retained  (in  percent) 


90% 


too-; 


Figure  9:  Percent  error  in  moment  at  mast  and  antenna  connection  plotted  versus 
the  percent  of  available  mast  modes  retained.  (Forcing  Frequency:  8.95  Hz) 


44 


Percent  error  in  moment  at  mast  Percent  error  in  shear  at  mast  and 

and  antenna  connection  antenna  connection 


Percent  error  in  shear  at  mad  and  antenna  connection  plotted  versus  the  percent  of 
available  antenna  modes  retained 


2  00% 

1  80%  * 

1  60%  • 

1  40%  - 
t  20% 

1  00%  ■ 

0  60%  - 
0  60  V;  - 
0  40%  ' 

0  20%  «  . 
0  00  %  -  - 
20% 


30% 


- •  -  CkAlG-BAMkH  iN 

-  -•  •  CRAki-CUANu 

a  MAC\1  \l 


------  ♦  -■=  . - 

40  Vo  50%  60%  70%  80% 

Number  of  antenna  modes  retained  (in  percent) 


90  % 


100  Vt 


Figure  10:  Percent  error  in  shear  at  mast  and  antenna  connection  plotted  versus 
the  percent  of  available  antenna  modes  retained.  (Forcing  Frequency:  8.95  Hz) 
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Figure  11:  Percent  error  in  moment  at  mast  and  antenna  connection  plotted  versus 
the  percent  of  available  antenna  modes  retained.  (Forcing  Frequency:  8.95  Hz) 
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igure  12:  Percent  error  in  shear  at  mast  and  antenna  connection  plotted  versus 
the  percent  of  available  mast  modes  retained.  (Forcing  Frequency:  219  Hz) 
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Figure  13:  Percent  error  in  moment  at  mast  and  antenna  connection  plotted  versus 
the  percent  of  available  mast  modes  retained.  (Forcing  Frequency:  219  Hz) 
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Figure  14:  Percent  error  in  shear  at  mast  and  antenna  connection  plotted  versus 
the  percent  of  available  antenna  modes  retained.  (Forcing  Frequency:  219  Hz) 
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Figure  15:  Percent  error  in  moment  at  mast  and  antenna  connection  plotted  versus 
the  percent  of  available  antenna  modes  retained.  (Forcing  Frequency:  219  Hz) 
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B.  BASE  EXCITATION  FROM  PRESCRIBED  DISPLACEMENT 


The  following  derivation  is  applicable  to  the  base  excitation  problem,  where  base 
displacements  (as  opposed  to  base  accelerations)  are  to  be  prescribed  As  in  the  base 
excitation  from  prescribed  acceleration,  the  derivation  starts  with  the  FEM  generated  mass 
and  stiffness  matrices  as  follows 


The  bottom  row  is  expanded  into  the  following  equation: 


+Kb^,X1.  +KBliXb  0 


From  Eq  (28).  the  base  acceleration  is  obtained  as  follows 

X|;  +Km.Xt,  +KBBX|4] 


(25) 


(28) 


(34) 


From  Eq  (25),  the  top  row  is  expanded  to  obtain 

Mou*o  "^MyBxB  +^oc>xc.  4"KOBxB  F0  (26) 

Substituting  Eq  (34)  into  Eq  (26)  and  simplifying,  the  equation  of  motion  for  the  internal 
coordinates  as  a  function  of  prescribed  base  displacement  is  obtained  as  follows 


[Moo  M0uMBl,M boJx0  *f[1x00  MobMbi,Kbo  Jxg  F0  [Kob  4"MobMbbK.bbJxb  (35) 
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Notice  that  the  matrix  that  pre-multiplies  the  acceleration  term  has  units  of  mass,  and  the 
matrix  that  pre-multiplies  the  displacement  term  has  units  of  stiffness  The  matrix  that 
pre-multiplies  the  base  displacement  term  has  units  of  stiffness  Therefore,  all  terms  of  this 
equation  of  motion  are  dimensionally  consistent  w  ith  units  of  force 

The  same  modal  decomposition  procedures  described  in  the  base  acceleration 
formulation  apply  to  the  base  displacement  formulation  Since  the  displacements  at  the 
base  are  prescribed,  the  base  coordinates  are  no  longer  degrees  of  freedom  As  in  the  base 
acceleration  formulation,  the  time  history  of  the  base  displacement  is  taken  to  be  simple 
harmonic,  and  is  represented  by 


{\11}={X,  }sin(a7)  (36) 

Numerical  convergence  examples  are  provided  in  the  following  section  Both  examples 
are  similar  to  the  examples  presented  in  the  prescribed  base  acceleration  problem 

1.  TIP  DEFLECTION  CALCULATION 

The  mast  and  antenna  system  was  subjected  to  a  base  excitation  where  the  time 
history  of  the  displacement  of  the  base  coordinates  was  prescribed  The  excitation  was 
performed  at  a  frequency  which  corresponded  approximately  to  mode  5  (47  26  Hz)  of  the 
mast/antenna  system  Since  the  system  was  modeled  without  damping,  a  frequency  which 
corresponded  exactly  to  a  natural  frequency  of  the  system  could  not  be  prescribed  The 
excitation  frequency  is  equivalent  to  the  mode  5  natural  frequency  to  within  2  decimal 
places  As  in  the  base  acceleration  problem,  the  frequency  was  arbitrarily  selected  The  tip 
deflection  of  the  antenna  was  calculated  and  percent  error  in  tip  deflection  was  plotted 
versus  percent  of  available  component  modes  retained  (see  figures  16  and  17)  The 
calculations  were  performed  twice  Mast  modes  were  truncated  in  the  first  calculation, 
and  antenna  modes  were  truncated  in  the  second  calculation  All  three  procedures  yielded 
accurate  results  in  both  the  mast  and  antenna  truncation  tests  When  the  mast/antenna 
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system  was  excited  at  mode  5,  the  Craig-Chang  and  MacNeal  procedures  provided  very 
similar  results  for  both  the  mast  truncation  and  antenna  truncation  As  was  determined  in 
the  tip  deflection  calculations  of  the  base  acceleration  problem,  the  Craig-Bampton 
procedure  yielded  the  best  results  using  fewer  modes  than  did  tlu  Craig-Chang  and 
MacNeal  procedures  in  both  truncation  tests  This  again  could  possibly  be  due  to  the 
combination  of  fixed  interface  normal  modes  and  static  constraint  modes  providing  a 
better  representation  of  prescribed  base  displacement  than  the  combination  of  free 
interface  normal  modes  and  residual  flexibility  modes  However,  all  three  provided 
accurate  results  and  at  a  cheaper  cost  than  the  FE  model 
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Figure  16:  Percent  error  in  antenna  tip  deflection  plotted  versus  the 
percent  of  available  mast  modes  retained  (Forcing  Frequency:  47.26  Hz) 
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Figure  17:  Percent  error  in  antenna  tip  deflection  plotted  versus  the 
percent  of  available  antenna  modes  retained.  (Forcing  Frequency:  47.26  Hz) 
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2.  MOMENT  AND  SHEAR  CALCULATION 

Percent  error  in  shear  and  moment  at  the  mast  and  antenna  connection  were 
plotted  versus  the  percent  of  av  ailable  component  modes  retained  (see  figures  18-21)  The 
base  excitation  was  conducted  again  at  47.26  Hz  Similar  results  were  obtained  in  the 
shear  and  moment  calculations  as  were  obtained  in  the  tip  deflection  calculations  The 
Craig-Bampton  procedure  provided  more  accurate  results  using  fewer  modes  than  the 
Craig-Chang  and  MacNeal  methods  The  convergence  rate  of  the  Craig-Chang  and 
MacNeal  methods  were  almost  identical  in  the  mast  mode  truncation  When  the  mast  was 
excited  at  mode  5,  all  three  yielded  excellent  results  in  the  antenna  mode  truncation 
analysis  However,  the  Craig-Bampton  procedure  converged  more  quickly  than  the  other 
methods  It  appears  that  the  combination  of  fixed  interface  normal  modes  and  static 
constraint  modes  lead  to  a  higher  rate  of  convergence  in  determining  tip  deflection  and 
antenna/mast  shear  and  moment  calculations  Although  all  three  procedures  initially  had  a 
higher  percentage  error  when  truncating  mast  modes,  than  when  truncating  antennae 
modes,  they  all  provide  as  accurate  if  not  more  accurate  results  than  the  antenna 
truncation  at  less  cost  in  terms  of  computations  It  is  also  noteworthy  to  compare  the 
computational  cost  in  retaining  mast  modes  versus  retaining  antenna  modes  in  predicting 
accurate  system  response  For  example,  in  the  shear  calculation  w'hen  the  mast  was 
excited  at  mode  5,  the  Craig-Chang  procedure  yielded  a  percentage  error  of  2  1%  while 
using  1,66  -lO6  FLOPS  during  the  mast  mode  truncation  test  During  the  antenna  mode 
truncation  test,  the  Craig-Chang  procedure  yielded  a  percentage  error  of  2.03°/o  using 
1.89  -106  FLOPS 

Since  the  impact  of  shock  waves  on  the  mast/antenna  system  are  typically  of  low 
frequency,  and  as  can  be  seen  from  the  results  of  the  prescribed  base  acceleration  and  base 
displacement  examples  provided,  it  is  recommended  that  the  Craig-Bampton  component 
mode  representation  be  used  to  synthesize  the  mast  and  antenna  system 
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Figure  18:  Percent  error  in  shear  at  mast  and  antenna  connection  plotted  versus 
the  percent  of  available  mast  modes  retained.  (Forcing  Frequency:  47.26  Hz) 
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Figure  19:  Percent  error  in  moment  at  mast  and  antenna  connection  versus 
the  percent  of  available  mast  modes  retained.  (Forcing  Frequency:  47.26  Hz) 
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Figure  20:  Percent  error  in  shear  at  mast  and  antenna  connection  plotted  versus 
the  percent  of  available  antenna  modes  retained.  (Forcing  Frequency:  47.26  Hz) 
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Figure  21:  Percent  error  in  moment  at  mast  and  antenna  connection  plotted  versus 
the  percent  of  available  antenna  modes  retained.  (Forcing  Frequency:  47.26  Hz) 
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VI.  PRELIMINARY  RESI  LTS  AND  CONCLUSIONS 


The  survivability  of  shipboard  combat  systems  equipment  is  paramount  to  the  warfare 
fighting  capability  of  the  ship  and  her  crew  Should  a  fire  control  radar,  or  any  vital 
topside  combat  systems  equipment  fail  as  a  result  of  an  induced  shock  wave,  the  ship's  w  ar 
fighting  capacity  would  be  crippled  However,  using  proven  structural  dynamics 
techniques,  the  design  engineer  can  design  the  mast/antenna  system  in  such  a  manner  as  to 
minimize  risk  of  failure 

As  seen  in  this  report,  the  mast  and  antenna  can  be  treated  as  separate  substructures, 
and  using  CMS,  can  be  assembled  as  a  masfantenna  system,  from  which  dynamic 
response  to  base  excitation  can  be  calculated  Treating  each  antenna  as  a  substructure, 
allow  s  the  cataloging  of  the  various  antennae  A  selected  antenna  can  be  "plugged"  into 
various  locations  along  the  mast  until  a  suitable  dynamic  response  is  obtained  As  has 
been  demonstrated  herein,  CMS  in  conjunction  with  FE  modeling  provides  rapid  and 
accurate  results  at  a  computational  cost  significantly  less  than  standard  FE  modeling  The 
mock  mast/antenna  system  used  in  this  report  consisted  of  only  17  elements  and  a  total  of 
51  degrees  of  freedom  Although  the  results  that  were  obtained  on  this  "small"  model 
were  accurate  and  required  small  compute  times,  the  same  benefits  demonstrated  herein 
can  be  expected  with  larger  models,  specifically  those  that  will  be  used  to  represent  real 
mast  and  antennae  systems 

All  three  methods  yield  results  that  are  accurate  and  more  computationally  efficient 
than  standard  FE  modeling  However,  from  the  results  obtained,  it  is  recommended  that 
the  Craig-Bampton  component  mode  representation  be  used  to  synthesize  the  mast  and 
antennae  system  As  mentioned  previously,  the  mast  is  excited  typically  by  low 
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frequencies  Since  the  Craig-Bampton  procedure  yielded  more  accurate  results  while  using 
fewer  component  modes  than  the  other  methods,  the  Craig-Bampton  procedure  is  the 
substructure  coupling  formulation  of  choice  due  to  good  accuracy  and  ease  of 
implementation  Additionally,  the  NASTRAN  superelement  scheme  contains  the  Craig- 
Bampton  component  mode  representation  as  a  solution  path  to  dynamic  analyses 
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